Quantum statistical model of nonlinear inverse bremsstrahlung absorption 

in strongly coupled plasmas 
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A new approach for the calculation of collisional inverse bremsstrahlung absorption of laser light 
in dense plasmas is presented. Quantum statistical formalism used allows avoiding ad hoc cutoffs 
that were necessary in classical approaches. Thus, the current method remains reliable for strong 
electron-ion interactions. In addition, both the dynamic, field dependent response and hard electron- 
ion collisions, are consistently incorporated. The latter were treated in an average manner as a 
stopping power that in turn was cast into a form of a friction force. Here, for the first time a link 
between the stopping power and the problem of collisional laser absorption is drawn. This allows 
the theories developed for the stopping power calculation, such as the quantum T-matrix approach, 
to be applied to the problem of collisional laser absorption. The new approach accommodates the 
low- and high-frequency limits explained in the text and is valid for arbitrary laser field intensities. 
A comparison with classical MD simulation is indicative of the validity of the new method in the 
wide parameter range tested. 

PACS numbers: 52.38.Dx, 52.38-r, 52.27.Gr 
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I. INTRODUCTION 

Understanding the laser-matter interaction in strongly 
coupled plasmas is crucial for the design of the con- 
temporary inertial confinement fusion (ICF) targets. In 
both the direct and indirect drive ignition schemes, the 
laser energy is deposited into the plasma of high-Z ele- 
ments such as Au. The conditions imposed by the hy- 
drodynamic instabilities on the spatial symmetry of the 
laser energy deposition, require the laser absorption to be 
carefully determined from the early stages of the laser- 
plasma interaction In the case of the fast ignition 
[2| it is important to model the laser absorption start- 
ing with the preformed plasma due to the unavoidable 
nanosecond prepulse. The critical plasma density cor- 
responding to the third harmonic of the Ny:Yag lasers, 



used in most inertial fusion designs is 
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In the limit of low temperatures (few cV) such plas- 
mas are characterised by strongly coupled and degen- 
erate electrons: T = (e 2 /fc B T e )(47rn e /3) 1 / 3 ~ 1 and 
n e (2,Trh 2 /m e kBT e ) 3 ~ 1 respectively. 

The dominant mechanism of radiation absorption for 
lasers with intensities typical for ICF is the inverse 
bremsstrahlung. In this case, the radiation is absorbed 
via collisions between the plasma particles usually de- 
scribed in terms of the electron-ion collision frequency 
Q. First calculation of the inverse bremsstrahlung in 
the high frequency limit for the lowest order in the laser 
field strength was reported by Dawson & Oberman (DO) 
Later Decker et al. Q have extended this result 
to arbitrary field strengths. A classical ballistic model 
was considered in Ref. [fj in order to study the frequency- 
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dependent electron-ion collisions in plasmas. 

These approaches arc formulated in the high-frequency 
limit, when the number of binary electron-ion collisions 
per laser cycle can be neglected. In this limit, the laser 
is coupled to the plasma via the induced polarisation 
current so that the electron-ion interaction has a collec- 
tive rather than a binary character. In the low-frequency 
limit strong scattering due to the binary collisions dom- 
inate the process of the absorption and the induced po- 
larisation becomes relatively small. At intcrmidiate fre- 
quencies the effects of binary and collective scattering 
have to be considered simultaneously. 

For a given laser frequency, plasma conditions can be 
such, that the rate of the binary collisions becomes com- 
parable with the laser period. Considering 0.35 /im lasers 
used in most applications, this occurs in a strongly cou- 
pled plasmas. The collisions in such plasmas have to 
be evaluated quantum mechanically. The first quantum 
treatment was reported in Ref. 0- The nonlinear ab- 
sorption was determined using the first Born approxi- 
mation in Rcfs. H @- A quantum approach to cal- 
culate laser absorption in strong fields was also devel- 
oped in Ref. [Tbl Semiclassical approach in the linear 
regime using a memory function kinetic formalism in- 
cluding lowest order quantum effects was developed in 
Ref. . A quantum statistical approach for dynamical 
conductivity in strongly coupled regime was developed 
in Refs. [lj, [HI]. A quantum Vlasov approach for ar- 
bitrary field strength, similar to the classical approach 
developed by Decker et al. [1[ was presented in Ref. 0- 
Rigorous kinetic approach to the inverse bremsstrahlung 
absorption in strongly coupled plasmas using nonequi- 
librium Green's function techniques was developed in a 
series of publications [H, [lr3l . The expressi on for the col- 
lision frequency derived by iBornath et all using the lat- 
ter approach is identical to the quantum Vlasov method 



2 



fl4| . In f act, it can be obtained from the formula ob- 
tained by iDecker et al.l by replacing the classical dielec- 
tric function by the quantum Lindhard dielectric function 
EH- This indicates, that both approaches make sim- 
ilar approximation, neglecting the effect of strong scat- 
tering by binary collisions, which makes them applicable 
in the high-frequency limit only. The effect of the bi- 
nary collisions was considered in a linear-response theory 
by using the Gould-DeWitt scheme in Ref. [H. Inverse 
bresstrahlung absorption for strongly coupled plasmas 
was also calculated using classical molecular dynamic 
(MD) simulations reported in Refs. [20|, [U [HJ . 

In this paper, we present a description of collisional 
absorption that bridges between the high- and low- 
frequency limits. It is shown that the interactions can 
be split into a weak collective interactions and hard col- 
lisions. The latter arc treated as the stopping power of 
ions in the electron fluid that can be cast into the form 
of a friction between the electron and ion fluids. Thus, 
for the first time, the stopping power formalism is ap- 
plied to the calculation of the collisional absorption, al- 
lowing one to use a well developed models of the stopping 
power (see, e.g., Refs. [ISHjEl]) m tne problem of laser 
absorption in plasmas. The description of the collective 
electron response can be kept almost unchanged from ear- 
lier approaches 0, [3] . Due to the use of full quantum 
mechanical formulation of the problem, no ad hoc cut- 
offs must be introduced and the theory stays reliable for 
strong electron-ion interactions and degenerate electrons. 
The few assumptions made are justified by the unprece- 
dente d ag reement with molecular dynamic (MD) simula- 
tions [20T[2ll. |22| up to very high coupling strengths. 

In the following chapters we develop a quantum formal- 
ism to describe the laser absorption in dense plasmas. In 
SecUJa set of Vlaso y-Poisson (VP) e quations similar to 
the one developed bv lKull and Plagnel is introduced. The 
major improvement in the present model is the inclusion 
of the hard-collisions and its treatment using the stop- 
ping power formalism. In Sec. IIIII the VP equations are 
solved in the Kramers-Hcnncbcrgcr (KH) frame that is 
calculated with respect to the effect of the hard-collisions 
on the electron fluid rest-frame. Results and discussion 
follow in Sec. IP71 



II. 



QUANTUM VLASOV EQUATION WITH AN 
EFFECTIVE FRICTION FORCE 



The motion of the electrons in a neutral plasma con- 
sisting of Ni ions and N e = ZNi electrons, where Z is 
the average charge state is described by a one electron 
statistical operator pi (t) whose evolution is governed by 
equation: 



Ho, p\ 



NiVi%p?o}, 



(1) 



The ions are considered to be located at fixed positions 
Xj, j = 1,2, . . . , Ni distributed according to the temper- 



ature dependant ion- ion pair correlation function ga(x). 
The Hamiltonian of the system is 



H 



K 

2m 



(2) 



where U ex t describes the externally applied field, p| e 2 and 
pl l 2 are the two particle electron-electron and electron- 
ion distribution functions respectively, and ~V\\ are 
the e-e and e-i interaction potentials. Next assuming that 
the two-particle density functions can be expressed as: 



~a3 
Pi, 2 



P?pt 



■912 



(3) 



The first term in the expansion of p^ f on the r.h.s is the 
Hartrec term, while all the higher order terms are con- 
tained in the the corresponding pair correlation function 

9 1,2- 



Using these definitions Eq. JT]) can be written as 



where 



H,pi 



^{NeV^ + NtV^}, (4) 



H = H - e$, (5) 
and the effective Hartrec potential $ is: 

- e$ = Ti{N e V^ 2 pF 2 + NiVftft) (6) 

In Eq. j4]) the hard, binary, collisions are grouped in the 
r.h.s. and weak collisions are treated in the framework of 
Hartrec approximation as a collective average potential 
contributed by the system species. This effective poten- 
tial is determined self-consistently from the charge den- 
sities of electrons and ions by the Poisson equation: 



Ni 

A $(x, t) = Ane{n e (x, t) - Z^^Sfa - xj 



(7) 



where n e (x,t) = N e (x\pi\x) is the electron density. The 
potential is calculated as a classical field ${x, t) which is 
used as an operator &(x,t) in Eq. (H}. Thus, assuming 
electrostatic interactions with a self consistent collective 
scattering potential §{x, t) and an externally applied po- 
tential $ ext = —x-E ext (t), due to a time dependent laser 
field E ext (t) = Eo sin(a;oi), the Hamiltonian can be writ- 
ten as 



Pi 



H = — e [$(aji,i) 
Am 



<5>ext{Xi,t)} 



(8) 



The heating rate of the plasma by the external laser 
field can be expressed using the effective electron-ion col- 
lision frequency 0] given by the following expression: 



Anul (j)-E _ 4mj2 



EE 



(9) 
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Here, the overline stands for averaging over one oscil- 
lation period and the angular brackets denote the ex- 
pectation values of the quantum operators. In the last 
equation j is the electric current density operator: j = 
— (en e /rn)pi. In order to determine the time dependence 
of the expectation values of the operators, one has to 
switch from the microscopic quantities given by the ki- 
netic equations to the average macroscopic quantities via 
the statistical operators. Thus, multiplying by pi and 
applying the trace to the both sides of Eq. ((4]) , the time 
change-rate of the expectation value of the momentum 
(pi) is obtained 



*h-( Pl ) 



Pi,H 



ihTrT? c pi, 



where the hard-collision integral 1^ is defined by: 



T? c ee ^{NiVi e ; 2 gf; 2 + N e V^gt%}/ih. 



(10) 



(11) 



The integral on the right-hand side contains hard- 
collisions only and is related to the hard-collisions con- 
tribution to the stopping power by 



dE 
d.r 



Tr{l^(t)pi). 



(12) 



In order to obtain an analytical solution of Eq. ([9]), it 
is useful to cast the effect of the stopping power as an 
average friction force between the electron and ion fluids 
R = —(dE/dx)/V, where V = ( P i)/m is the average 
ensemble velocity, and restrict the solution to small par- 
ticle velocities, such that V < Uthi where v t h = kT e /m 
is the electron thermal velocity, since in this case the 
friction coefficient is velocity independent [U Hi| . The 
dynamics of screening is of minor importance in the low 
velocity range, justifying the use of the stopping power 
data calculated assuming statically screened Coulomb in- 
teractions [HlUll. In this limit, the rate equation for the 
average momentum becomes 

-(pi) = -e[E sin (w„t) + (E)] - ^ c (pi), (13) 

where the hard collision frequency ^hc is defined as v\ xc = 
R/m and the polarisation field is (E) = (d$>/dx). 

Field amplitudes for which the assumption of velocity 
independent friction coefficient R is valid are restricted 
by the condition Vmax < Vth- Neglecting the contribution 
of the polarisation field (E) in Eq. (|13p . the maximum 
average velocity can be estimated by 



V m 



1 + "he 



1 + 



(14) 



he 



where vq = eEo/muJo is the free electron quiver velocity, 
and v\ ic = v\ ic /uj a is the normalised hard-collision fre- 
quency. Therefore, in the high-frequency limit: i\ c <C 1 
the fields amplitudes are restricted by vq < fth! and in 



the low-frequency limit: Phc 3> 1 the condition becomes 
Vq < Phc v th- Thus, because of the strong damping of 
the velocity in the low- frequency limit, the region where 
R = const is applicable can be extended to stronger 
fields. 

Furthermore, the use of the stopping power formalism 
imposes additional limitation with respect to the laser 
frequency loq. The stopping power treatment assumes 
that the typical collision time is greater than the period of 
plasma oscillation characterised by the plasma frequency 
Lu p ee Aire 2 n e /m. Therefore, the analytical solution to 
be obtained is restricted to the laser frequency range of 

Multiplying both sides of Eq. (fl~3|) by —en e /m the cur- 
rent balance equation is obtained 

^ = ^ [E sin (w t) + (E)] - v hc (j) (15) 

Note, that in the low- frequency limit ^> 1, the polar- 
isation field (E) vanishes and Eq. (fT"3f produces the well 
known Drudc formula for the low-frequency conductivity: 



4it<t d = v - . 

^hc - «w 

A formal solution of Eq. (|15|) is 



— Ul 'Y Eq 

ti)(t) = ~ ~r [ cos (uot) - v hc sin (uj t)} 

47T UJq 

j (E)(r)e^-^dr 



(16) 



(17) 



where 7 = 1/(1 + v£ c )- One can see that the current 
consists of two terms, namely the polarisation current 
represented by the second term on the r.h.s., and the 
free current, represented by the first term. The latter 
includes a contribution from the hard collisions that are 
in phase with the laser field. The collision rate is obtained 
by substituting Eq. ([17]) into Eq. ©: 



2 ^o 1 

: 7 I/ hc + -^2" I E sm{u a t) 



\E)(t)< 



(18) 



Here, the first term is due to the the strong collisions. 
It has the same form as the real part of the Drude con- 
ductivity. The second term is due to the polarisation 
current. 

To complete the derivation one needs to calculate the 
polarisation field (E) to be used in the last equation. To 
do so the system of Eq. (01 and Eq. ([7]) has to be solved. 
The general solution of this problem is notoriously diffi- 
cult however, it can be greatly simplified if one assumes 
that the same average friction force — Vhc(Pi) is acting on 
all the electrons irrespective of their direction, position 
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and velocity. This assumption is justified in the limit of 
u th n > wnen the ions directed motion as a particle 
beam characterised by a single velocity with respect to 
the electrons can be considered. In such case, Eq. can 
be cast into the form of the quantum Vlasov equation: 



lH di Pl 



H e s,pi 



0, 



(19) 



where the effective Hamiltonian 7J e ff is: 

H eS = |i- - e ($ + $ ex t ) + Vhc(pi)x! (20) 

Eq. (fl"9]l produces the same current balance equation (fT~5|) 
as Eq. {4}. The set of two equations: Eq. {7} an Eq. HH) 
form a closed set of Vlasov-Poisson (VP) equations with 
the hard collision determined using the standard methods 
applied for the stopping power calculations. 



in Ref. [Tj and skip some parts of the derivation elabo- 
rated in Ref. [Til while underlining the differences obtained 
in the end result. 

In the KH frame, the scatterring potential is treated 
as a perturbation H^' = — e$ c ff of the equilibrium set 
up by the Hamiltonian = p 2 /2m. Correspondingly, 
the density operator is set to be p = p^ ' + p^ l \ where 
p^ ) is the unperturbed stationary ensemble. 

The Fourier-Laplace transform of the effective poten- 
tial Eq. (|22bp is found to be: 



D(k,uj) 



(23) 



where £i(fc,cj) is the Fourier-Laplace transform of the 
ion potential in the KH frame. The Fourier transform 
assumes the form, 



E 1 -(fc,t) = S il0 (fc) x e lfc '^ 



(24) 



III. SOLUTION OF QUANTUM VP 
EQUATIONS 

In this chapter we shall obtain the solution of the sys- 
tem of VP equations introduced above. The laser and the 
friction act both as an effective external force since they 
depend on x only. Therefore, in the absence of the scat- 
tering field $ the electrons would perform a quiver mo- 
tion. It is useful to transform to Kramers-Henneberger 
(KH) reference frame - the rest frame of the electron 
fluid, since one can assume that in this frame the elec- 
trons are close to the equilibrium unperturbed state. The 
transformation to KH frame is given by y — x — and 



where 



£ = -je [sin (u> t) + v hc cos (u t)] 



(21) 



where e = —eEo/muj 2 .. In this frame the VP equations 
become: 



p 

2m 



e$ KH ,p KH 



Ni 



= 



(22a) 



A $ KH (a;, t) = 47re(n e (y, t) - Hv + (22b) 



Equations (|2Tj) and ([22)1 form the basis for the calcu- 
lation of the laser energy absorption by the electron-ion 
scattering. Since in the following analysis we shall carry 
out the calculation mainly in the KH frame, we shall 
omit the KH superscript. The solution procedure of this 
set of equations in identical to the solution procedure of 
the VP equations (4a) and (4b) discussed in Ref. 0- The 
only difference is in the definition of the KH frame. Here, 
the effect of the hard collisions in the form of the friction 
force is included in the determination of the KH frame, 
whereas in Ref. [U the electrons are freely oscillating in 
the external field. Therefore, we adopt the notations used 



(25) 



is a static part and 

oo 



is a dynamic phase factor due to the quiver motion. Its 
Fourier coefficients are Bessel functions of the first kind 
J n depending on the parameter z = k • e. Applying the 
Laplace transformation to Eq. (|2"4")l yields 



Ei(fe,w) = £*,(,(*) 2 



(-l)"(-») ffl J B (7z) J m (7^hc z) 
—i(uj + (n + m)oJo) 



(26) 

The electric potential generated by the ions in the KH 
frame is screened by the dielectric function D(k,uj) 



D(k,uj) = 1 + 4irx(k,uj) 



(27) 



which is known as the Lindhard dielectric function [LJ 
HH and the response function x(k,u>) is given by 



X(k,u>) = ^ / d 3 u 



e 



f(u + Hk/2m) - f(u - Hk/2m) 



h{uj — k ■ u) 
The perturbed charge density 

r.(fc,w) = -&N eP W(k,u) 



(28) 



is related to the effective potential and the response func- 
tion by 



T e {k,u) = -k 2 X (k,Lu)<Z>(k,Lu) 



(29) 
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Finally, wc apply the inverse Laplace transform to 
T e (fe, u) and uj) and obtain from the poles of (k, uj) 
at the frequencies uj = — (n + m)wo the asymptotic result, 



*(M) 



(fe)e 



■i(n+m)ct;o^ 



(30a) 



(30b) 



where the coefficients of these series are given by 

" (-l)"(-i) m J„( 7 z) J m ( 7 P hc z) 



$n,m(fc) = £ i>0 (fc)^ 



L>(fc, -(n + m)wo) 



(31a) 



T n , m {k) = -k 2 x(k, -(n + m)w )$n, m (fc) (31b) 

These results allow us to obtain the expectation value 
of the electric field used in the calculation of the colli- 
sion frequency in Eq. (|18[) . In the momentum representa- 
tion E(k, t) = —ik§(k,t) and using the transformations 
(A8b), (A15) from Ref. Qj yields 



(E(k,t)) 



1 



d 3 k 



-ik$*(k,t)T e (k,t) (32) 



Substituting this expression into Eq. (fTS)) and using 
Eqs. (|30|) we get the collision frequency 



UJ V 



w p A"/ fc2 D(k,(n + m)u ) 

n— — oo 

oo 

^ « V m+S (7^2) (n - s)(l + i Vhc) Jn-s (7 Z ) 
s— — oo 

- i (7^c z)J n -B-i iiz) 5ii(fe) 

where, Sa(k) is the ion-ion structure factor, and 
a = (2^ 2 )- 1 (w /w P )(^e 2 /m«2) 



(33) 



(34) 



Clearly, the known limiting cases can be readily retrieved. 
Decker's result (Eq. 20 in Ref. |]| follows from Eq. (33]) 
in the weak coupling limit i>i ic —> and non-degenerate 
plasmas. In this case D(k,uj) becomes the classical di- 
electric function and the integral must be truncated at 
kmax to avoid the divergence at small impact parameters 
(see Ref. [2?] for the discussion of different cuttofs). Here, 
all integrals can be performed to infinity and no ad hoc 
cutoffs must be introduced as a result of the quantum 
mechanical treatmen. The first term in line 1 dominates 
for small laser frequencies ?hc ~^ 00 giving a Drude-likc 
expression. 



IV. RESULTS AND DISCUSSION 

In the previous sections we have presented a quantum 
mechanical formulation for the problem of laser absorp- 
tion in dense plasmas. The for malism is similar to that 
developed by iKull and Plagnel . but inherently includes 
the hard-collisions absent from this and other approaches 
[S EE EH- This was achieved in three steps. Firstly, 
the the two-particle density function was split into: (i) 
the first order Hartree term, representing the weak in- 
teractions; (ii) the higher orders, representing the hard- 
collisions, which were collected into the pair correlation 
function (e.g. Eq. (J3J). Secondly, the hard-collisions were 
cast into the form of the friction force using the stopping- 
power integral on the right-hand side of the momentum 
rate equation (JTUJ) . At last, the collision frequency Eq. @ 
was determined by finding the first order perturbation of 
the equilibrium density distribution as a result of the 
electron-ion scattering. The equilibrium electron distri- 
bution is set up in the electron KH rest frame determined 
by the external potential and the friction force due to the 
hard-collisions. 

The use of the stopping power formalism restricts this 
approach by demanding that: (i) the interaction time is 
larger than or at least of the same order of magnitude 
as the typical plasma oscillation period, setting ujq < c^ p ; 
(ii) the field strength is sufficiently high to consider the 
ion motion relative to the electrons as that of a directed 
beam characterised by a single velocity, i.e. V max 3> v 1 ™. 
From the other hand, the analytical solution of Eq. © 
expressed in Eq. (|33[) is only valid in the low-velocity 
limit V max < «th, where the friction coefficient R is con- 
stant. The value of the maximum relative velocity V max 
depends both on the laser field strength and Vh c accord- 
ing to Eq. HID. 

The end result in Eq. (f33)) formally resembles the 
Gould-DeWitt ansatz [28|| due to the splitting of the col- 
lision frequency into a sum of two contributions resulting 
from the strong and weak interactions. However, here no 
ad hoc assumption was made and the splitting to the hard 
and weak collision contributions in Eq. (J33J was obtained 
as a result of the discussed solution process. Thus, the 
present approach might also hint on the region of appli- 
cability of the Gould-dcWitt scheme when used in other 
models. 

The hard electron-ion collisions are incorporated in 
Eq. ([33]) via a friction force related to the stopping power 
of the ions in an electron gas that in turn sets up a 
more general KH frame instead of the freely oscillating 
one adopted in other works 0, [f| [lj, [l5[. Many mod- 
els have been developed for the stopping power [231 ] . few 
include hard collisions. Within quantum statistical the- 
ory, they can be described by a T-matrix approach based 
on the quantum Boltzmann equation. The related cross 
sections are calculated from numerical solutions of the 
Schrodingcr equation [24|. The full stopping power can 
be then determined applying the Gould-DeWitt scheme 
[28| or by velocity-dependent screening length [25|, [26| . 
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Pfalzner & Gibbon [20]: MD 
Hilse et al. [21]: MD 
Morozov et al. [22]: MD 
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Bornath et al. [16] 
Decker et al. [5] 
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FIG. 1: Collision frequency u ei versus coupling parameter F 
for a hydrogen plasma with fixed temperature / density and 
a laser field with ojo/uj p = 1 (upper panel), uiq/uj p = 3 (lower 
panel) and «o/fth = 0.2. Solid line: Eq. (|33[) ; punctured lines: 
contributions of hard collisions and polarisation to Eq. (|33|) . 
The classical results of Decker et al. H) was calculated with 
an integral cut-off at k rnax = rn e v^ h /Ze 2 . 



According to the conditions of Fig. (Q]), the first approach 
is sufficient and has been used to generate the data pre- 
sented here. That is, the hard collision term has been 
calculated as the stopping power using full cross sections 
(T-matrix approach) minus the one in static first Born 
approximation. 

Next, the numerical solutions of Eq. ([55)) are com- 
pared with the earlier approaches. We start with a 
commonly accepted be nchmark case first considered by 
iDawson and Obermanl and later used in works by dif- 
ferent authors, ft is presented in Fig. [1] where our re- 
sults ( e.g . Eq. (|33|) ) arc compared with other theories 
IE 03, Q El and simulation data [H HH, H2 . Yet an- 
other reason for choosing this particular regime is that it 
presents the hardest parameter set for comparison since 
it is lying inbetween the known limiting cases. As ex- 
pected, all theories agree for weakly coupled plasmas, 
but large deviations occur for strong coupling. The clas- 
sical description is clearly not a pplicable here as demon- 
strated by comparison with the iDecker et alj result (also 
see Ref. [27| for extended discussion). For a coupling 
strength of T > 1, the quantum theories of Refs. [lj, [l5| 
also start to disagree with the simulation data. Such 
discrepancy can be traced back to the neglect of hard 
collisions. Here this shortcoming is overcome. As a re- 
sult we find an excellent a gree ment with data from MD 
simulations by Iflilse et "all and M orozov et all Despite 



FIG. 2: Real part of the conductivity as a function of the laser 
frequency for a hydrogen plasma n e = 10 22 cm -3 , vo/v t h = 
0.2. 



the inconsistency between these two data sets, our result 
reproduces the main features observed in both simula- 
tions, especially the sharp increase of v ei for T > 1 rela- 
tive to other theories [HI, HH and the change in the sign 
of the slope. Plotting both contributions of Eq. (|33|) sep- 
arately reveals that the hard-collisions term dominates 
for high coupling strengths and defines the shape of the 
curve in this parameter area. From this, the observed 
change in the slope sign is due to the turnover in the 
hard-collision contribution. The latter occurs because 
the static Born overtakes the T-matrix contribution at 
high coupling strengths. 

Note, that different parameter sets are tested in the up- 
per and lower frames of Fig. [1] the data running as a func- 
tion of the increasing density in the upper and decreasing 
temperature in the lower frame, and good agreement be- 
tween the MD results and our approach is obtained for 
this wide data range. However, the contribution of hard- 
collisions is rather small for the parameters presented in 
the upper frame, and therefore the advantages of the cur- 
rent approach are less pronounced in that case. 

Degeneracy might obscure the comparison at low tem- 
peratures or high densities. However, degeneracy is nei- 
ther included in the MD simulations nor in our calcula- 
tion of the hard collision term which is based on solutions 
of two-particle Schrodinger equation. We therefore com- 
pare our data to the MD simulations on a similar level 
of approximation. 

The multi-dimensional parameter space should also be 
examined along the direction of the laser frequency. This 
comparison is demonstrated in Fig. [2] for fixed electron 



7 



density and the ratio of vo/vth = 0.2, where the conduc- 
tivity is plotted as a function of the laser frequency. Here 
again we show the contributions of Eq. |33|) separately. 
As expected, the contribution of the dynamic, polarisa- 
tion term vanishes at low frequencies luq <C to p , and the 
conductivity is dominated here by the Drude-like term 
due to the hard-collisions. At the intermediate frequency 
range u>o ~ u p , both terms contribute equally, with the 
dynamic term overtaking for lower values of T. As dis- 
cussed earlier, the present approach is only valid for the 
laser frequencies that are of the same order of magnitude 
as the plasma frequency and lower, therefore we do not 
extend the comparison to the high laser frequencies. The 
breakdown of the present approach at high frequencies is 
already exhibited in Eq. (|33|) . As follows from the latter, 
at the high frequency limit the hard-collision dominate 

the absorption, moreover v e i 0> P > v\ lc = const, which 
is obviously sensless. 

In conclusion, a quantum mechanical approach for the 
calculation of collisional absorption of laser light in dense 
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